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ABSTRACT 

We discuss the dependence of the pulsation frequencies of the axial quasi- normal 
modes of a nonrotating neutron star upon the equation of state describing the star 
interior. The continued fraction method has been used to compute the complex fre- 
quencies for a set of equations of state based on different physical assumptions and 
spanning a wide range of stiffness. The numerical results show that the detection of 
axial gravitational waves would allow to discriminate between the models underly- 
ing the different equation of states, thus providing relevant information on both the 
structure of neutron star matter and the nature of the hadronic interactions. 

Key words: gravitational waves - equation of state - quasi-normal modes of compact 
stars. 



1 INTRODUCTION 

The search for gravitational waves has received a big burst 
during the last decade. In addition to the already oper- 
ating resonant antennas (EXPLORER, AURIGA, NAU- 
TILUS, ALLEGRO, NIOBE), ground-based interferometers 
(VIRGO, LIGO, GEO600, TAMA) will be soon operating, 
and a future space-based interferometer, LISA, will extend 
the exploration to low frequencies, unaccessible to other de- 
tectors. Thus, the detection of gravitational waves is close to 
become reality. Among the astrophysical processes that are 
associated with emission of gravitational waves, the gravi- 
tational collapse and the coalescence of compact bodies are 
likely to be the most efficient sources. Both processes leave 
behind a compact body, either a neutron star or a black hole, 
which is expected to radiate the mechanical energy residual 
of its violent birth through gravitational waves. These waves 
will be emitted at frequencies and with damping times char- 
acteristic of the quasi-normal modes of the source. It is in- 
teresting to note that these modes may be excited also in 
other astrophysical processes that can be observed. For in- 
stance, stellar oscillations are thought to be at the origin 
of the drifting subpulses and micropulses detected in some 
radio sources, and of the quasi-periodic variability seen in 
some X-ray burst sources and in a number of bright X-ray 
sources (McDermott, Van Horn & Hansen 1988). 

The complex frequencies of the quasi-normal modes 
carry information on the internal structure of the emitting 
source. For black holes, it has been shown that they depend 
exclusively on the parameters that identify the spacetime ge- 



ometry: the mass, charge and angular momentum. For stars, 
the situation is far less simple, since the quasi-normal mode 
eigenfrequencies depend on the equation of state (EOS) pre- 
vailing in the interior, on which not much is known. It is 
therefore interesting to compute these frequencies for differ- 
ent EOS's proposed to describe matter at supernuclear den- 
sities, and explore the possibility of extracting from them 
information on the internal structure of the star. 

The quasi-normal mode frequencies can be computed by 
studying the source-free, adiabatic perturbations of an equi- 
librium configuration with an assigned EOS, and by solving 
the linearized Einstein equations, coupled with the equations 
of hydrodynamics, with suitably posed boundary conditions. 

If the unperturbed star is assumed to be static and non 
rotating, it is convenient to expand all perturbed tensors 
in tensorial spherical harmonics, and since these harmonics 
have a different behaviour under the angular transforma- 
tion 6 TV — 6, ip ~* TV + ip, the separated equations 
split in two decoupled sets: the polar or even, belonging to 
the parity ( — 1)^, and the axial or odd, belonging to the 
parity (—1)'^^^'. The polar equations are the relativistic 
generalization of the tidal perturbations of newtonian the- 
ory, and couple the perturbations of the gravitational field 
with the perturbations of the fluid composing the star. As 
in newtonian theory, the classification of the polar modes is 
based on the behaviour of the perturbed fluid according to 
the restoring force that is prevailing: the g-modes, or gravity 
modes, when the force is due to the eulerian change in the 
density, the p-modes, when it is due to a change in pressure. 
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and the f-mode, that is the generaUzation of the only pos- 
sible mode of oscillation of an incompressible homogeneous 
sphere (Chandrasekhar 1964). This classification scheme was 
introduced by Cowling (Cowling 1941). In addition to the g, 
f, p modes, that exist also in newtonian theory, in general 
relativity there exists a new family of modes that are es- 
sentially spacetime modes, since the corresponding motion 
of the fluid is negligible (Kokkotas & Schutz 1992). They 
are named w-modes, and are characterized by frequencies 
typically higher than those of the g, f, p modes, and much 
smaller damping times, i.e. these modes are highly damped. 

The frequency of the f-mode was determined by Lind- 
blom and Detweiler for several EOS's (Lindblom & Detweiler 
1983), and more recently this study has been extended to in- 
clude a few modern equations of state suggested for neutron 
stars, for which also the frequencies of the first p-mode and 
of the fundamental polar w-mode have been computed. In 
addition, it has been shown that it is possible to infer empir- 
ical relations between the mode frequency and the macro- 
scopic parameters of the star: the mass and the radius (An- 
dersson & Kokkotas 1998). 

Unlike the polar perturbations, the axial perturbations 
are not coupled with fiuid motion. They do not have a new- 
tonian counterpart, and their radial evolution is described 
by a Schrodinger-like equation with a potential barrier that 
depends on the distribution of energy and pressure in the 
interior of the star in the equilibrium configuration (Chan- 
drasekhar & Ferrari 1991a). Thus, the EOS of the fluid has 
the sole role of determining the shape of the potential inside 
the star. 

The axial quasi-normal modes divide in two classes: 
the w-modes, highly damped and with properties similar to 
the polar w-modes, and the slowly damped s-modes (Chan- 
drasekhar & Ferrari 1991b). These modes appear if the star 
is extremely compact, so that the potential well in the inte- 
rior becomes deep enough to allow for the existence of one 
or more quasi-stationary states, i.e. of quasi normal modes. 

In this paper we shall calculate the frequencies of the 
axial quasi-normal modes for the following EOS's; whenever 
possible, we follow the classification scheme introduced by 
Arnett & Bowers (1977): 

• Pandharipande, model A. The star is described in 
terms of pure neutron matter, the dynamics of which 
is dictated by a nonrelativistic hamiltonian containing 
a semi-phenomenological interaction potential (Pandhari- 
pande 1971a). 

• Pandharipande, model B. A generalization of model A 
including the appearance of protons, electrons and muons in 
/3-equilibrium, as well as of heavier baryons (hyperons and 
nucleon resonances) , at sufficiently high densities (Pandhari- 
pande 1971b). 

• Wiringa, Fiks and Fabrocini, model WFF. Neutron star 
matter is treated as a mixture of neutrons, protons, elec- 
trons and muons in /9-equilibrium. The nuclear hamilto- 
nian includes two- and three-body potentials. With respect 
to models A and B, the ground state energy of neutron star 
matter is computed using a more sophisticated and accurate 
many-body technique (Wiringa, Fiks & Fabrocini 1988) . 

• Akmal, Pandharipande and Ravenhall, model APRl. 
Similar to model WFF, but uses state of the art parametriza- 



tions of both the two- and three-body potentials (Akmal, 
Pandhripande & Ravenhall 1998). 

• Akmal, Pandharipande and Ravenhall, model APR2. 
Includes selected relativistic corrections to model APRl 
(Akmal et al. 1998). 

• Pandaripande & Smith 1975, model L. Neutrons are 
assumed to interact through exchange of uo — , p— and 
(T— mesons. While the exchange of heavy particles {uj and p) 
is described in terms of nonrelativistic potentials, the effect 
of the (7— meson is taken into account using relativistic field 
theory and the mean-field approximation (Pandharipande 
& Smith 1975). The physical assumptions underlying this 
hybrid approach are somewhat different from those of the 
previous models. 

A useful way of classifying EOS's is through their stiffness, 
which can be quantified in terms of the speed of sound Vs : 
stiffer EOS's correspond to higher sound speeds. The stiff- 
ness is also a measure of compressibility: stiffer EOS's cor- 
respond to more incompressible matter. 

The EOS's considered in this paper can then be ordered 
according to increasing stiffness as follows: B, A, APR2, 
APRl, WFF, L. 

Though the s-modes have been shown to exist for homo- 
geneous stars Qwith high enough compactness, the EOS's we 
have considered do not lead to potential wells deep enough 
to allow for the appearance of any s-mode. Thus we have 
focussed on the imprint that different EOS's leave on the 
axial w-modes. In particular, we have computed the pulsa- 
tion frequency and the damping time of the fundamental 
w-mode for the aforementioned EOS's, by using the contin- 
ued fraction method, first developed by Leaver to find the 
quasi-normal mode frequencies of black holes (Leaver 1985, 
1986), and which we have suitably modified to make it ap- 
plicable to the present context. 

We shall show that, as for the polar w-modes, it is possi- 
ble to derive empirical relations that allow to determine the 
mass and radius of the star. However, the most remarkable 
result is that the range within which the pulsation frequen- 
cies vary is peculiar of each EOS, much more than in the 
polar case. 

The plan of the paper is the following. In section 2 we 
shall write the equations governing the axial perturbations 
of a non rotating star. In section 3 the EOS's will be shortly 
described and compared. In section 4 the continued fraction 
method will be outlined and section 5 will be devoted to a 
discussion of the results. 



2 A SCHRODINGER EQUATION FOR THE 
AXIAL PERTURBATIONS 

The equations for the axial perturbations of a non rotat- 
ing star can be considerably simplified by introducing, after 
separating the variables, a function Zi(r), constructed 
from the radial part of the axial metric components, where 
I is the harmonic index (Chandrasekhar & Ferrari 1991a). 



* It should be mentioned that e = const, e being the energy 
density, is the stifTest EOS, since Vs = oo. 
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Ze{r) satisfies the following equation 

Vt{r)]Zi = 0, 



Ve{r) 



+ \UJ 



(1) 



-[;(/ + l)r + r-^(e-p)-6m(r)], 



where 



-V + IJ.2 



dr, and e and p are the energy 



density and the pressure of the perfect fluid composing the 
star. The function u{r) appearing in the potential barrier, 
can be found by integrating the equation 

= (2) 

and imposing the boundary condition that, at the surface 
of the star, the metric reduces to the Schwarzschild metric, 
i.e. e^"^^^ = 1 - where M = m{R) is the mass of 

the star. Outside the star, where e and p vanish, eq. (|l|) 
reduces to well known Regge- Wheeler equation 

Vi{r)]Zt^Q, (; 



dr-f 



Vtir) = —[l{l + l)r- 



6A/], 



and 



2v 



^ 2A; 
r 



Thus the axial perturbations of a star are fully described 
by a Schrodinger-like equation with a potential barrier that 
depends on the distribution of energy-density and pressure 
inside the star in its equilibrium configuration. 

The quasi-normal modes are defined to be solutions of 
eq. which satisfy the condition of being regular at r — 0, 
and that behave as a purely outgoing wave at radial infinity, 
i.e. 

Ze ~ g-'-^^* j.^ _> oo. (4) 



3 MODELS OF NEUTRON STAR MATTER 
EQUATION OF STATE 

In this work, we have used six different models of neutron 
star matter EOS at supernuclear density (e > 2.8 • 10^* 
g/cm^). To obtain the EOS over the whole density range 
relevant to neutron stars, models A, B , WFF and L have 
been matched to the EOS at subnuclear densities obtained 
by Feynman, Metropolis & Teller (1949), Baym, Pethick 
& Sutherland (1971) and Baym, Bethe & Pethick (1971), 
whereas models APRl and APR2 have been supplemented 
with the recent results by Pethick, Ravenhall & Lorenz 
(1995). However, the details of the EOS at subnuclear den- 
sity have negligible infiuence on both the neutron star prop- 
erties and the calculated pulsation frequencies. The equa- 
tions of state we consider are plotted in fig. 1 over the entire 
range of density, whereas in fig. 2 the different high density 
behaviours are compared. 

In fig. ^ we plot the mass radius relation for the above 
mentioned EOS's, in the range where we have checked that 
the stellar models are stable against radial perturbations. 

The EOS's denoted by A (Pandharipande 1971a), B 
(Pandharipande 1971b), L (Pandhripande et al. 1975) and 
WFF (Wiringa et al. 1988) have been widely used in the past 
to study equilibrium configurations of neutron stars (Arnett 
et al. 1977; Pandharipande, Pines & Smith 1976; Wiringa 
et al. 1988) and, more recently, to compute their pulsation 
frequencies associated with emission of gravitational waves 




log (p) 

Figure 1. The logarithm of pressure (dynes-cm^) versus the log- 
arithm of the energy density (g-cm~^) for the EOS's considered 
in this paper over the whole density range. In the region of sub- 
nuclear densities the curves corresponding to different EOS's sit 
on top of each other. 
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Figure 2. A blow-up of the region of supernuclear densities of 
figure ^ The units are the same as in the previous figure. It can 
be seen that, in the region relevant to neutron star properties, 
the EOS's can be ordered according to increasing stiffness as ex- 
plained in the text: B, A, WFF, APR2, APRl, L. 



(Andersson & Kokkotas 1998). On the other hand, the mod- 
els denoted APRl and APR2 (Akmal et al. 1998) have been 
developed within the last two years and have never before 
been employed to study nonstatic neutron star properties. 

With the only exception of model L, the EOS's consid- 
ered have been obtained using nonrelativistic nuclear many- 
body theory. Within this approach, nuclear matter is treated 
as a system of pointlike protons and neutrons, whose dynam- 
ics is described by the nonrelativistic hamiltonian: 
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i j>i k>j>i 

where m and Pi denote the nucleon mass and momen- 
tum, respectively, whereas Vij and Vijk describe two- 
and three-nucleon interactions. The two-nucleon potential is 
obtained from a fit to the properties of both the bound and 
scattering states of the two-nucleon system (Wiringa, Stoks 
& Schiavilla 1995) , while the three-body term Vijk has to 
be included in order to account for the binding energies of 
the three-nucleon bound states (Pudliner et al. 1995) . 

Obtaining the EOS of nuclear matter at zero tempera- 
ture amounts to evaluate its ground state energy per particle 
Eo as a function of baryon density, which in turn requires 
the solution of the many-body Schrodinger equation with 
the hamiltonian of eq.(|^. The EOS's considered in this work 
have been obtained using the Rayleigh-Ritz variational prin- 
ciple and the formalism of cluster expansions, which allows 
to write Eq as the sum of contributions arising from sub- 
systems containing an increasing number of particles. In the 
pioneering works of Pandharipande (1971a, 1971b) the cal- 
culations were carried out at lowest order in the cluster ex- 
pansion, i.e. including only two-body cluster contributions. 
The more recent WFF, APRl and APR2 EOS's have been 
obtained using the Fermi-Hyper-Netted-Chain (FHNC) in- 
tegral equations (Fantoni & Rosati 1975, Pandharipande 
& Wiringa 1979) to take into account the contributions of 
the relevant terms in the cluster expansion to all orders. 
It has to be pointed out that the FHNC equations allow 
for a very precise calculation of the ground state expecta- 
tion value of the nuclear hamiltonian. Moreover, within this 
computational scheme it is possible to use highly realistic 
trial wave functions, whose structure reflects the complexity 
of the nuclear interactions. In this respect, the results ob- 
tained by Wiringa et al. (1988) and by Akmal et al. (1998) 
are more accurate than those obtained by Pandharipande 
(1971a, 1971b). 

As far as the nuclear hamiltonian is concerned, it is im- 
portant to realize that the two-body potential models devel- 
oped over the past few years (Wiringa et al. 1995), as the Ar- 
gonne nig potential used in model APRl and APR2, rep- 
resent a significant improvement over the previously avail- 
able ones. They are fitted to a larger and more complete 
database of nucleon-nucleon scattering data and explicitly 
include the effect of charge-symmetry-breaking terms in the 
scattering amplitude. Results by Akmal et al. (1998) show 
that using a two-body potential of the last generation leads 
to qualitative changes in the density dependence of nuclear 
and neutron matter energy, suggesting the possibility of a 
transition from the standard uniform phase to a spin-isospin 
ordered phase associated with condensation of neutral pious 
(Tamiya & Tamagaki 1988, Benhar 1981). 

We will now briefly summarize the main features of the 
model EOS's used in this work. In the case of EOS A (Pand- 
haripande 1971a) neutrons are allowed to interact through 
a two-body potential, the three-body potential is neglected 
and the cluster expansion is truncated at the two-body level. 
The hamiltonian has been constructed using an early model 
of Vij, developed by Reid in the late sixties (Reid 1968). 

A more realistic description of neutron star matter re- 
quires inclusion of the effect of inverse /3-decay, which even- 
tually leads to the appearance of a mixture of neutrons, pro- 
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Figure 3. The mass-radius relation is plotted for the EOS's con- 
sidered in this paper, in the range of stability. M is expressed in 
solar mass units, and _R in km. 



tons, electrons and possibly muons in /3-equilibrium. Models 
B, WFF, APRl and APR2 are aU based on this picture. As 
already mentioned, the calculations by Wiringa et al. (1988) 
and Akmal et al. (1998) have been carried out using the 
FHNC equations. Among the different EOS's discussed in 
the first reference, we have used the one corresponding to the 
so called Urbana v\4 two-body potential, supplemented 
by a purely phenomenological three-nucleon interaction. 

From the computational point of view, the EOS of 
model APRl (Akmal et al. 1998) has been obtained us- 
ing basically the same approach as Wiringa et al. (1988). 
The main difference between models APRl and WFF comes 
from the nuclear hamiltonians. Model APRl uses the new 
Argonne v\g, potential and a three-body interaction, referred 
to as Urbana IX, featuring an attractive component, due to 
pion-exchange processes with excitation of a A resonance 
in the intermediate state, and a phenomenological repulsive 
part at shorter range. 

The relevance of relativistic corrections to the standard 
nuclear many-body picture has been studied in model APR2 
(Akmal et al. 1998) by including boost corrections to the 
nucleon-nucleon interaction of order (v/c), v being the ve- 
locity of the center of mass of the interacting pair. In order to 
fit the observed properties of the three-nucleon systems af- 
ter the inclusion of relativistic corrections, the change in the 
two-body part of the interaction requires the three-body po- 
tential to be modified accordingly. The resulting three-body 
interaction is referred to as UIX*. 

The main problem associated with models A, B, WFF, 
APRl and APR2 is that they lead to a violation of causality, 
since the predicted speed of sound exceeds the speed of light 
at very high density. However, this pathology only occurs at 
densities much greater than those relevant to the study of 
neutron star structure. 

The above problem can be avoided treating neutron star 
matter within relativistic field theory, with a lagrangian de- 
scribing the nuclear interaction in terms of one-boson ex- 
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change processes. The coupling constants are chosen in such 
a way as to reproduce selected nuclear properties, and the 
corresponding field equations are solved in the mean-field 
approximation. Following the pioneering work of Walecka 
(1974) this approach has been used by many authors. Model 
L (Pandharipande et al. 1975) includes the exchange of two 
vector mesons {uj and p) and one scalar meson (a). It can 
be regarded as a somewhat hybrid model, since only the a- 
exchange interaction is treated within the mean-field ap- 
proximation, whereas the exchange of the heavier vector 
mesons is described by static potentials. 

The EOS's employed in this work can be classified ac- 
cording to their stiffness, providing a measure of the in- 
compressibility of neutron star matter. Comparison between 
models A and B shows that the appearance of heavy baryons 
leads to a significant softening of the EOS at high density. 
This behaviour can be easily understood noting that, due 
to their low concentrations, the heavy baryons carry a ki- 
netic energy much smaller than that carried by the neutrons. 
Hence, replacing neutrons with heavy baryons produces a 
sizeable decrease of the pressure. 

In the case of neutron matter, either pure or in /3- 
equilibrium, the stiffness of the EOS is mainly determined 
by the nuclear hamiltonian. The EOS's of models A, WFF, 
APRl and APR2, in which the nuclear hamiltonian is con- 
trained by nucleon-nucleon scattering data, are rather close 
to each other. On the other hand, using a hamiltonian that 
does not reproduce nucleon-nucleon data, as in model L, 
may result in a significantly difference in stiffness. 

It is interesting to note that the comparison between 
models APRl and APR2 suggests that the inclusion of rel- 
ativistic corrections leads to a softening of the EOS. Hence, 
the fact that model L provides the stiffest EOS should be as- 
cribed to the different dynamical assumptions, rather than 
to its relativistic nature. 



4 THE CONTINUED FRACTION METHOD 

As stated in section ^ the axial quasi-normal modes of a 
pulsating star are solutions of eq.(0) satisfying the boundary 
conditions imposed by physical requirements: Zi should 
be regular at the origin, and have the behaviour of a purely 
outgoing wave at infinity. To compute the mode frequencies 
one has to find the (complex) values, uj = ujq -\- iun , for 
which these boundary conditions are satisfied. 

The numerical determination of the quasi-normal mode 
frequencies is non-trivial, especially for modes with large 
imaginary parts (strongly damped modes). The reason is 
simple to understand. Solutions of eq.(0) representing out- 
going and ingoing waves at infinity have the asymptotic be- 
haviour 



and 



-r* /t 



(6) 



as r« ^ oo, where r — l/cui is the damping time. 
Therefore, identifying by numerical integration the purely 
outgoing solutions (that is, those solutions for which Z}" 
is zero) becomes increasingly difficult as the damping of the 
mode increases. 

The same problem occurs also in the case of quasi- 
normal modes of black holes, and was solved by Leaver 
(Leaver 1985, 1986). He found a continued fraction relation 



that can be regarded as an implicit equation which identifies 
the quasi-normal frequencies, thus circumventing the need 
to perform an integration out to big values of r«. This 
method was subsequently reformulated and applied to the 
polar oscillations of a star (Leins, NoUert & Soffel, 1993). 
Here we show that it can be used also in the case of axial 
perturbations. 

In this section we will use dimensionless geometrical 
units, c = G = 2M = 1. In these units the Regge- Wheeler 
equation, which describes the perturbed spacetime outside 
the star, becomes: 



d^Zt 



+ 



drl 
and 

r* = r + ln(r — 1) 



0, 



(7) 



(8) 



We shall now write the solution of the Regge- Wheeler equa- 
tion in a power-series form as follows. Defining v = 1 — ^, 
where r = a is some point outside the star, and introducing 
a function <l>i{v), related to Zi{r) by: 

Z,{r) = (r - l)-"-e-"->,(«) = x{r)Uv), (9) 

one finds that (f)t satisfies the differential equation: 

{co + ClV + C2v'^ + C-j,V^)(t)l^^^ + (10) 

4-(do + dxv + d2V^)(j)i^^ + (eo -I- eiv)4>t = 0. 

The constants depend only on u, £ and a through the 
relations: 



Co = 1 

a 

do — —2iuja 



ci = - - 2 
a 

3-2 
a 



C2 = 1 - - 

a 



eo 



C3 

d2 



ei 



a 
3 
a 

3 

a 



Let us now perform a power expansion of (l>i{v): 



(11) 



By substituting this expression in eq.(^o|), the expansion co- 
efficients a„ are found to satisfy a four-term recurrence 
relation of the form: 



a„an+i + Pnan -\- 7„a„_i -I- 5„a„_2 
where: 



0, 



n>2, (12) 
(13) 



a„ = n(n + l)co 
Pn — {n— l)nci + ndo 
'Jn = {n- 2)(n - l)c2 + (n - l)di + eo 
Sn = {n- 3)(n - 2)c3 + {n - 2)^2 + ei 

The coefficients ao and ai can simply be determined by im- 
posing the continuity of Zi and Zi^r in r = a, since from 
eq.(p|) it follows that: 



ao 



ai 



Ze{a) 
Xia) 
a 



(a) + -—^Zi{a) 



(14) 
(15) 
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The values of Zi{a) and Ze ,.[a) can be obtained by 
numerically integrating eq. ( 0) in the interior of the star and 
continuing the solution outside, up to r = a, by integrating 
the Regge- Wheeler equation. The remaining coefficients can 
then be determined by recursion from eq.(|l2). 

To apply the continued fraction technique one has to 
deal with three-term recurrence relations. Leaver has shown 
that the four-term recurrence relation ( p^ can be reduced 
to a three-term relation by a gaussian elimination step 
(Leaver, 1990). In other words, with the positions (n = 0, 1): 



and for n > 2: 



1-71 — <-tn , 



7n = In 



7n-l 

eq.jl^) reduces to: 



7n-l 
<5n =0, 



(16) 



(17) 



(18) 



The elimination step is not as trivial as it may seem, be- 
cause in the process one of the three independent solutions 
to eq.(]l^) is lost. It can be shown (Leins et al. 1993) that 
this solution is not relevant for our purposes. 

We now turn to investigating the asyrnptotic behaviour 
in the expansion ([ll[). Let us make 



of the coefficients 
the ansatz: 



lim 



h 



k 

= i + -r72 +- + ■ 



(19) 

-3/2 



Dividing eq. ( |l^ by n^a„, keeping terms up to ~ n 
and equating to zero the various terms in the expansion in 
powers of n~^^'^ we find the relations: 



Co + Ci + C2 + C3 = 0, 
2co + Ci - C3 = 0, 

= 2iija, 
3 



(20) 



+ iLo{a + 1). 



The first two of these equations are identities. Substituting 
the second pair of equations in ^ we get: 



lim a„ = n 

n — 'oo 



-3/4+iu(a + l)g±2V2iuari 



(21) 



According to a definition given by Gautschi, the solution of 
eq. ( p^ ) corresponding to the plus sign in eq. (|2^) is said 
dominant, whereas that corresponding to the minus sign is 
said minimal (Gautschi, 1967). If we select the minimal so- 
lution we see that the expansion ( pT| ) is absolutely and uni- 
formly convergent outside the star, provided that we choose 
a such that 



a/2 < R < a 



and 



a > 2, 



(22) 



and that, according to equation (|9|), the solution to equation 
(^) behaves as a pure outgoing wave at infinity, i.e. it is the 
wavefunction of a quasi-normal mode. Thus, the key point is 
to identify the minimal solutions of eq. (|ig|). According to a 
theorem due to Pincherle (Gautschi 1967), if equation ( [l9| ) 
has a minimal solution then the following continued fraction 
relation holds: 



Oi 

ao 



-71 Q!l72 a?273 



(23) 



Pi- 132- Pz- 

where the continued fraction on the RHS is convergent and 
completely determined since the coefficients q„, and 
7„, defined in eqs. ( p^) , (p^ are known functions of u) . 
Moreover, from eqs. (Q and (|l|) it is apparent that the 
dependence on the stellar model is all contained in the ratio 
ai/ao. 

Defining now: 

A ai „ 
po = — , ao = -1, 
ao 

eq.(E3) can be recast in the form: 



(24) 



= h{oj) = l3o ■ 



ao7i ai72 0273 
/32- Pi- ' 



(25) 



Using the inversion properties of continued fractions (Wall 
1948), the latter equation can be inverted n times to yield: 



= /„(a;) = /3„ 



Qn-l7n an-27n-l ao7l 



I3n-1- 

an7n+l Qn+l7n+2 Qn+27n+3 
Pn + 1— Pn + 2— Pn+3 — 



/3n-2- 

for 



Po 



(26) 



These n conditions are analytically equivalent to eq.(p5[). 
Anyway, since the functions fn{i^) have different con- 
vergence properties, each of them is best suited to find the 
quasi-normal frequencies in a given region of the complex 
ui plane. This is the main reason behind the power and 
flexibility of the continued fraction technique. 

To find the mode frequencies we have adopted the fol- 
lowing numerical procedure. We compute the real and imag- 
inary parts of fn{(^) for a given inversion index n (say 
n = 0) on a suitably chosen grid of {u!o,(^i) values. We 
plot the curves along which the two functions are zero and 
look for the intersections of those curves: the quasi-normal 
frequencies as the points where the lines cross. For each as- 
signed value of the inversion index n, besides the physi- 
cal frequencies this method singles out some spurious roots, 
which, however, can be ruled out easily, since for different 
values of n they disappear. 

The method has been tested by verifying and extend- 
ing results on constant density models of stars previously ob- 
tained by Kokkotas (1995), who found the QNM frequencies 
through a direct integration of eq.(|l]). Our results are in ex- 
cellent agreement (up to the fifth significant digit) with those 
of Kokkotas. Furthermore, exploiting the inversion proper- 
ties of continued fraction, we were able to find an additional 
set of modes, the so called w//-modes (Leins et al. 1993), 
with still larger imaginary part. The real and imaginary part 
of the axial modes of a homogeneous star with R/M — 2.4, 
found with the continued fraction method, are plotted in fig. 
^ and compared with the quasi-normal mode frequencies of 
a Schwarzschild black hole. 



5 NUMERICAL RESULTS 

The continued fraction method has been applied to com- 
pute the frequency and the damping time of the first axial 
w-mode, Vwq and r^g, respectively, for models of stars with 
the EOS's discussed in section 3. In table 1 we summarize 
the characteristics of the considered stellar models. For each 
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Figure 4. The real and the imaginary part of the frequency of the 
w-modes (square) and of the w/j-modes for a homogeneous star 
with R/M = 2.4, are compared with the quasi-normal mode 
frequencies of a black hole. The diamond indicates the s-modc of 
the star. 



EOS, we choose the upper value of the central density (col- 
umn 2) as that of the last radially stable configuration. The 
lower value is chosen as in (Andersson & Kokkotas 1998) for 
the EOS's A, B, WFF and L, whereas for APRl and APR2 
- that have been used to compute the oscillation modes of 
a star in this paper for the first time - it is chosen to give a 
value of M/R comparable with that of the other models. 
In column 3 we give the ratio M/R corresponding to the 
selected endpoints of pc, the radius and mass of the star 
are given in columns 4 and 5, respectively, and the values of 
and r™, are in columns 6 and 7. 



In fig. H we plot 



as a function of the star compact- 



ness, M/R, for the EOS's given in table 1 and within the 
corresponding central density range. The damping time, Twq , 
is similarly plotted in fig. ^. Since the pressure and density 
inside the star are obtained by logarithmically interpolating 
the pressure-density tables, for values of M/R close to the 
stability limit their reconstructed radial behaviour is less ac- 
curate. Hence the corresponding values of u^^ and r^^ are 
less reliable. For this reason, we have omitted from table 1 
and fig. I and I the marginally stable configurations for the 
WFF equation of state. 

From table 1 and fig. ^ we see that, apart from the 
equations of state WFF, APR2 and APRl, for which the 
ranges within which varies partially overlaps, the 

frequencies of the lower w-mode for the remaining EOS's, i.e. 
B, A, WFF and L range within intervals that are separated. 
This means that if an axial gravitational wave emitted by a 
star at a frequency of, say, 10 kHz could be detected, then 
we would be able to identify the equation of state prevailing 
in the star interior as the EOS B (among those we have 
considered), even without knowing the mass and the radius! 

In order to check whether a similar identification can be 
done by detecting a polar gravitational wave, it is interesting 
to compare the frequencies of the fundamental axial w-mode 
with those of the fundamental polar w-mode. The latter have 



Table 1. For each EOS, we give the minimum and maximum val- 
ues of the central density (g-cm""^) chosen as explained in section 
5 (column 2). In correspondence of these endpoints we tabulate: 
the ratio M/R in geometrical units (column 3), the radius (km) 
and mass {Mq) of the star (columns 4 and 5, respectively), the 
values of (kHz) and t^q (ps) of the lowest axial w-mode 
(columns 6 and 7). 




APRl 



A — AA L 



0,15 



0,20 0,25 
M/R 



0,30 



Figure 5. The frequency of the first axial w-mode is plotted as 
a function of the compactness of the star. 



been computed by Andersson and Kokkotas (Andersson & 
Kokkotas 1998) for several equations of state, including EOS 
A, B, WFF and L, which we have also considered. Using 
their results, in fig. ^ we plot the frequencies of the first axial 
and polar w-mode as a function of the star compactness for 
these EOS's, and in table 2 we tabulate the same quantities 
given in table 1, and the values of the polar u^^ and t^q 
corresponding to the central density range they consider. 
From table 2 and fig. ^ we see that the detection of a polar 
gravitational wave at a frequency of, say, 10 kHz, would 
not allow to discriminate between the equations of state A, 
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Figure 6. The frequency of the first polar (dashed Une) and 
ajcial (continuous Une) w-modes are plotted as a function of the 
compactness of the star for the EOS's A, B, WFF, L. 




0,15 



0,20 0,25 

M/R 



0,30 



Figure 7. The damping time of the first w-mode is plotted as a 
function of the compactness of the star. 
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Table 2. The central density, radius, mass, and the frequency 
and the damping time of the polar w-mode are tabulated, as 
in table 1 (same units), for four equations of state considered 
both by Andersson and Kokkotas and in the present paper. The 
data are extracted from the tables Al, A2, A9 and AlO given in 
(Andersson & Kokkotas 1998). 



give information on this stellar parameter. To this purpose, 
we have fitted the values of v^n^ and t^,^ for the six 
EOS considered in this paper, respectively with a linear and 
a quadratic function of the compactness M/R. We find 
the following robust empirical relations 

1 



R 



M ; 



(27) 



where is expressed in kHz, and R in km, with 

a = -155.45 ± 3.23, (28) 
b= 121.69 ±0.73, 

and ~ 1.84. The damping time is well approximated by 
the relation 



"0 



1 

M 



2 Af 
+ b—+c 



where 



R 

is in ms, and 



solar mass units, with 

a= -1672.74 ±63.39 
b = 482.57 ± 27.40 
c = 37.37 ±2.80. 

,2 



(29) 

Al is the mass expressed in 
(30) 



and X = 2.12. 



WFF and L, and similar considerations hold for different 
values of the wave frequency. 

Our results show that the frequencies of the axial grav- 
itational waves emitted by neutron stars oscillating in their 
w-modes are mainly driven by the stiffness of the EOS, 
which is in turn dictated by the physical assumptions un- 
derlying the model describing the star interior. As a con- 
sequence, models WFF, APRl and APR2, based on very 
similar assumptions (/3-stable matter with hadronic interac- 
tions constrained by nucleon-nucleon scattering data) yield 
frequency ranges that significantly overlap with each other. 

Although the data available on radiopulsars allow to es- 
timate the mass of the observed neutron stars, little is known 
on their radius. It is therefore useful to provide methods that 



6 CONCLUDING REMARKS 

Since the axial perturbations of non rotating stars are de- 
scribed by a Schrodinger equation with a potential barrier 
shaped by the manner in which pressure and density are dis- 
tributed inside the star, the frequencies of the axial quasi- 
normal modes are expected to depend on the equation of 
state. This was shown to happen also for the polar w-modes 
(Andersson & Kokkotas 1998) that are coupled to negligi- 
ble fluid motion. In particular, both the axial and the po- 
lar w-modes appear to depend essentially on the stiffness 
of the equation of state, and this is further manifested by 
the possibility of finding a robust fit of both the frequency 
and the damping time of these modes as a function of the 
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compactness of the star. However, for each selected EOS the 
frequency of the polar w-modes is a rather steeply decreasing 
functions of M / R, whereas in the case of the axial ones 
the dependence of u^^ on the star compactness turns out 
to be weak. This different behaviour is clearly displayed in 
fig. ^, where the frequency of polar and the axial wo-mode 
are compared. As a consequence of this behaviour, the axial 
w-modes give a more direct and explicit information, com- 
pared to the polar ones, on which equation of state prevails 
inside the star, regardless of its mass and radius. 

We emphasize that our results show that EOS's based 
on different assumptions correspond to non overlapping fre- 
quency ranges. Hence, the detection of axial gravitational 
waves would allow to further constrain the existing models, 
with regard to both the composition of neutron star matter 
and the description of the hadronic interactions. The real 
problem is whether we will ever be able to detect an axial 
gravitational wave impinging on earth from a ringing star. 
Only detailed simulations of astrophysical processes, like the 
gravitational collapse or the coalescence of compact bodies, 
will tell us to what extent can the w-modes be excited and 
be significant from the point of view of detection. If this will 
happen to be the case, since these modes have frequencies 
higher than those detectable by current experiments, new 
high frequency detectors will need to be planned. 
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